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摘 要 : 探讨 获 散 (87) 变 化 特征 及 其 影响 因素 对 区 域 水 资源 有 效 利用 至 关 重要 。 本 文 基于 BEPS(Boreal Ecosys- 


tem Productivity Simulator) fK AN ,结合 遥感 数据 和 和 气象 数据 模拟 了 内 蒙古 半 干 旱 区 1981 一 2018 年 57 变 化 特征 并 量 
化 了 其 影响 因子 的 贡献 。 研 究 发 现 ,1981 一 2018 年 内 蒙古 半 干 旱地 区 ET 以 1.75 mma 的 速度 呈 波 动 上 升 趋势 (P< 
0.05) ;同时 ,ET 变化 存在 年 代 际 差异 ,在 1997 年 发 生 明 显 突变 ,1998 一 2018 年 有 7 以 1.70 mm*a' 的 速度 增加 (P< 
0.05)。 大 气 水 分 气 缺 (VPD) 和 LAI 在 1997 年 后 每 年 分 别 以 0.002 hPa 和 0.01 的 速度 显著 增加 (P<0.05 ) ,弥补 了 其 他 
因子 的 不 利 影响 ,导致 了 ET 显著 增加 。 通 径 分 析 和 回归 分 析 结 果 表 明 ,ET 变化 主要 受到 VPD 的 直接 作用 驱动 , 通 


径 分 析 的 决定 系数 达 0.95, 由 VPD 主导 17 变化 的 区 域 


24.83%~90.46%。 同 时 ,VPD 也 是 耕地 、 林 地、 草地 ,城市 币 


能 够 解释 45% 以 上 的 57 变 化 。 
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ET 是 地 表 水 分 转移 到 大 气 的 过 程 , 是 连接 地 表 
水 和 能 量 平衡 的 关键 纽带 "。 由 于 气候 变化 和 人 为 
因素 影响 ,全 球 &7 整 体 旦 上 升 趋势 ;ZT7 变 化 将 导致 
水 循环 加 速 ” ,从 而 影响 区 域 水 资源 的 供需 情况 。 
因此 ,准确 模拟 和 量化 区 域 尺 度 ET 及 其 影响 因素 有 
助 于 进一步 理解 气候 变化 对 陆地 水 循环 的 影响 。 

在 过 去 的 几 十 年 中 ,已 经 对 7 的 变化 和 相关 的 
驱动 力 进 行 了 一 些 研究 ”; 但 由 于 影响 因子 变化 的 
空间 异 质 性 和 区 域 气候 背景 的 差异 ,目前 对 引起 ET 
变化 的 影响 因子 依然 没有 定论 ”。 近 年 来 ,大 气 水 
分 亏 缺 (VPD) 和 植被 作为 两 种 重要 的 影响 因子 受到 
WIET ZKE., Zeng Em T hist ETXE 
被 绿化 响应 的 最 新 文献 和 数据 集 , 并 指出 绿化 对 全 
球 陆 地 87 增加 的 贡献 达 50% 以 上 5 。 气 候 变 暖 导 
致 饱和 水 汽 压 增强 ,海洋 输送 水 汽 减 少 导致 大 气 实 
际 水 汽 压 降低 ,大 气 干旱 加 剧 ; 在 干旱 半 干 旱 区 ， 
持续 升温 将 加 剧 VPD 对 有 7 的 影响 呈 。 而 在 缺 水 、 
低产 的 草原 地 区 ,ET 对 这 种 气候 变化 的 影响 最 为 敏 
感 "。 内 蒙古 半 干 旱 区 以 草地 为 主要 植被 类 型 , 生 
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占 内 蒙古 半 干 旱 区 面积 的 93.56% ,能 够 解释 ET 变化 的 
1 裸 地 5 种 土地 利用 类 型 上 7 变化 的 主要 影响 因素 ,平均 


态 脆弱 ,对 气候 变化 和 人 为 干扰 敏感 ;水 资源 时 空 
分 布 不 均 且 亏 缺 严重 ,是 全 球 增 温 速 度 最 迅速 的 区 
域 之 一 ”。 目 前 ,对 陆地 储 水 量 等 的 分 析 均 显示 出 
内 蒙古 半 干 旱 区 存在 干旱 加 剧 的 趋势 ;因此 ,在 
当前 气候 背景 下 ,明确 内 蒙古 半 干 旱 区 的 57 变 化 及 
其 驱动 因子 以 应 对 水 循环 加 速 可 能 造成 的 极端 事 
件 , 对 区 域 水 资源 科学 调度 和 风险 管理 提供 决策 
依据 。 

本 文 以 长 期 通 量 和 气象 数据 为 基础 ,结合 
BEPS (Boreal Ecosystem Productivity Simulaton) 模 型 ， 
探究 1981 一 2018 年 内 蒙古 半 干 旱 区 ET 的 时 空 变化 
特征 ,揭示 其 对 气候 和 植被 变化 的 啊 应 ,以 期 量化 
ET 的 变化 特征 及 其 主要 影响 因子 ,对 区 域 水 资源 有 
效 开发 利用 和 生态 环境 保护 具有 重要 的 科学 意义 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 
依据 国家 气象 局 编 绘 的 中 国 气 候 区 划 图 为 标 
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准 ,选取 位 于 中 国 西北 部 的 内 蒙古 半 干 旱地 区 作为 
研究 区 域 (107.66°~121.90°E,37.62°~46.83°N)。 研 
究 区 为 典型 的 农 牧 交错 地 区 ,生态 环境 脆弱 ,是 对 
气候 变化 响应 最 为 敏感 的 地 区 之 一 。 草 原 占 总 面 
积 的 62% 以 上 ' ,主要 分 布 于 内 蒙古 高 平原 的 北 
部 ;9% 的 森林 主要 分 布 于 内 蒙古 高 平原 的 南部 和 大 
兴安 岭南 部 地 区 ;耕地 占 15% 左 右 ,分 布 于 内 蒙古 
高 平原 的 南部 和 松 辽 平原 西南 部 (图 1)。 
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图 1 研究 区 位 置 及 土地 利用 类 型 分 布 
Fig. 1 Study area location land use and cover type 
distribution in 2018 


1.0 数据 来 源 

本 文选 用 来 自 于 国家 气象 科学 数据 中 心 (http:// 
data.cma.cn) 的 1981 一 2018 年 内 蒙古 192 个 站 点 的 
气象 数据 ,数据 均 经 过 严格 的 质量 控制 和 标准 化 订 
正 处 理 , 主 要 包含 逐日 平均 温度 (TA,%C ) 降水 
(PRE,mm*a )\ 日 照 时 数 (SD,h)、 风 速 (WD ms) 
和 相对 湿度 (RH,%)。 大 气 水 分 亏 缺 (VPD ,hPa) 数 
据 集 由 上 述 气 象 数据 计算 获得 。 土 壤 湿 度数 据 来 
自 欧 洲 中 期 天 气 预 报 中 心 (https://cds.climate.coper- 
nicus.eu) 的 土壤 湿度 网 格 数据 集 ,该 数据 集 采 用 在 
2-5 cm 深度 的 表层 土壤 层 中 的 体积 土壤 水 分 
(SMV ,mm ) 来 表示 土壤 湿度 。 叶 面积 指数 (LAI1) 
数据 来 自 资源 环境 科学 与 数据 中 心 的 GLOBMAP 
LAI V3 数据 集 (https://www.resdc.cn/) ,本 研究 中 以 
年 均 LAI 代 表 当 年 植被 的 生长 状态 “”。 上 述 数据 格 
式 均 处 理 为 WGS1984 投 影 和 1 km 空间 分 辩 率 。 

1981 一 2018 年 1 km 空间 分 辩 率 的 87 数 据 集 是 
基于 日 步 长 BEPS 模 型 模拟 获得 。BEPS 模 型 是 在 叶 
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片 生理 生态 模型 FOREST-BGC 的 基础 上 建立 的 ,是 
一 种 基于 遥感 的 生态 机 理 过 程 模 型 ( 双 叶 双 源 模 
型 )。 模 型 中 7 计算 考虑 了 冠 层 截留 降水 蒸发 、 冠 
层 蒸 腾 、 土 壤 表 面 蒸 发 和 地 表 积 雪 升华 4 部 分 ;综合 
考虑 了 环境 (土壤 水 分 、 气 象 要 素 等 ) .植被 结构 和 
植物 生理 等 因素 对 ET 的 影响 。BEPS 被 证 实 能 够 准 
确 的 估算 ET 并 揭示 气候 变化 和 植被 动态 变化 对 7 
的 影响 机 制 ,因而 得 到 了 广泛 应 用 "*。 

ET 验证 数据 来 自 Asia Flux 观测 通 量 数据 库 
(http://www.asiaflux.net) ,选取 2006 一 2013 年 多 伦 和 
锡林浩特 两 个 典型 的 半 干 旱 草地 站 的 逐 月 通 量 数 
据 ,站 点 位 置 如 图 1 所 示 。 该 组 通 量 数据 的 能 量 平 
衡 比率 达 91% 以 上 ,符合 能 量 不 闭合 的 要 求 ,数据 
质量 良好 ”1。 

1.3 研究 方法 

1.3.1 ET 变化 特征 分 析 ”本文 从 突变 检验 和 趋势 
分 析 两 个 角度 来 描述 长 时 间 序 列 内 蒙古 半 干 旱 区 
Z7 变 化 特征 。 首 先 , 采 用 累积 距 平 法 (CA) 和 
STARS(Sequentialt-Test Analysis of Regime Shifts ) 两 
种 不 同 的 突变 检测 方法 检验 内 蒙古 半 干 旱地 区 
1981—2018 年 7 数据 突变 情况 ,两 种 方法 交叉 验 
证 可 以 有 效 的 提高 突变 检验 精度 。 选 择 两 种 突变 
含 验 中 相同 突变 年 份 作为 突变 点 ,对 有 7 变化 情况 进 
行 分 段 研究 。 累 积 距 平 法 (CA) 是 一 种 能 够 通过 累 
积 曲线 直观 地 判断 数据 演变 趋势 和 突变 检验 的 非 
线性 统计 方法 ”。STARS 是 由 Rodionov 提出 的 一 
种 基于 顺序 :检验 检测 序列 平均 波动 水 平和 方差 变 
化 的 方法 ,能 够 实时 监测 系统 从 一 种 相对 稳定 的 状 
态 迅 速 转变 为 另 一 种 状态 过 程 的 信号 ,该 方法 现 已 
广泛 应 用 于 气候 生态 系统 突变 检验 研究 ”。 突 
变 可 通过 格局 转变 指数 (Regime Shift Index, RSI) 
PS TE 


j*m x. 
RSI, = 3. m=0, 1, 1-1 (1) 
7 £e, 
X; = 一 天 (2) 
x o T Ep + diff (3) 


diff -1 zc (4) 


式 中 :1 为 设置 的 序列 x 的 截止 长 度 , 类 似 于 低 通 滤 
波 中 的 截止 点 ;假设 序列 x 中 数据 到 第 1 个 截止 长 度 
1 计 为 格局 及 ,平均 值 x 作为 RI 的 度量 ; o NA 
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止 长 度 ! 序 列 x 的 平均 方差 ;为 在 给 定 概率 水 平 下 
RA 21-2 自由 度 对 应 的 :分布 值 。 按 照 若 RSI 值 为 
负 ,在 1 年 不 存在 突变 ,重新 计算 z, 寻找 下 一 个 起 
始点 对 应 j 值 ;RSI 值 为 正 值 ,格局 在 j 年 发 生 突变 。 
采用 Theil Sen 和 Mann-Kendall 方 法 进行 £7 趋 
势 分 析 , 这 种 方法 是 一 种 非 参 数 方法 ,通常 被 用 来 
进行 水 文 -气象 时 间 序 列 趋势 分 析 关 区; 其 中 Theil 
Sen 方 法 用 于 计算 内 蒙古 半 干 旱 区 87 的 变化 率 ， 
Mann- Kendall 方法 用 于 检验 趋势 的 显著 性 。Theil 
Sen 方 法 中 ,ET 的 变化 率 B 计算 方法 如 下 : 
i=j 
式 中 : x 和 2 IS UR EFA x AARE AE B 20 , 
序列 呈 上 升 趋势 ; 若 B<0 ,序列 呈现 下 降 趋 势 。 
1.3.2 蒸 散 (ET) 影 响 因子 分 析 ET 变 化 是 由 太阳 
辐射 .空气 温度 .空气 湿度 .风速 ,水 供应 以 及 植被 
的 生理 特性 等 因素 共同 驱动 。 由 于 影响 因子 的 复 
杂 性 和 相关 性 ,首先 通过 偏 相 关 分 析 和 通 径 分 析 的 
方法 ,在 温度 .降水 .相对 湿度 .风速 太阳 辐射 E 
壤 湿 度 和 LAI 中 选择 合适 的 影响 因子 作为 回归 分 析 
的 自 变 量 。 
通 径 分 析 能 够 在 多 元 线性 回归 基础 上 将 自 变 
量 和 因 变 量 的 相关 关系 分 解 为 直接 影响 和 间接 影 
响 ,以 直接 通 径 系数 和 间接 通 径 系 数 表征 自 变 量 对 
因 变 量 的 影响 ,具体 计算 方法 见 文献 [25]。 本 研究 
中 , 因 变 量 为 87, kT 数据 通过 了 Shapiro-Wilk Test 
正 态 性 检验 ,服从 正 态 分 布 ,可 以 进行 通 径 分 析 。 
决定 系数 绝对 值 越 大 ,表明 自 变 量 对 因 变 量 影响 强 
度 越 强 ; 自 变量 的 决定 系数 为 正 时 ,对 因 变 量 起 到 
促进 作用 ;反之 为 抑制 作用 *。 


B= Median[ ) Vi»j (5) 


模拟 ET/mm 


y-1.0031x-0.20 
R1-0.81, P«0.05 


观测 ET/mm 


在 减少 影响 因子 间 相 关 性 影响 的 基础 上 ,基于 
逐步 多 元 回归 分 析 来 探究 7 变化 的 主要 影响 因 
子 。 逐 步 多 元 回归 分 析 方 法 被 用 来 量化 每 种 环境 
子 对 ET 变化 的 相对 贡献 ”。 首 先 ,对 原始 数据 和 
环境 变量 进行 归 一 化 : 


x' PE i min (6) 


式 中 , x' 是 归 一 化 变量 ; x, 是 序列 x 中 第 ;个 变量 ; 
xus 和 x 分 别 是 序列 x 中 最 小 值 和 最 大 值 。 通 过 
逐步 多 元 回归 分 析 建 立 模 型 描述 归 一 化 后 ET 和 影 
响 因 子 间 的 关系 ,其 中 7 作为 因 变 量 Y, 影 响 因子 
作为 自 变量 X。E7T 影 响 因 子 归 因 计 算 如 下 : 
Y-aX,*bX,*cX,* --- (7) 


g- lel 
加 DT 
式 中 :ab 和 cc 是 回归 系数 ,9 是 对 了 变化 的 相对 


2 结果 与 分 析 


2.1 ET 数据 集 验 证 

为 了 验证 BEPS 模 型 对 研究 区 域 87 的 模拟 能 
力 , 选 取 2006 一 2013 年 锡 林 格 勒 和 多 伦 两 个 站 点 逐 
月 ET 观测 数据 与 BEPS 模 型 模拟 的 ET 数据 集 进行 
了 对 比 验 证 (图 2)。 与 观测 的 ET 相 比 ,两 个 站 点 上 
BEPS 模 型 模拟 的 ET 结果 均 能 够 较 好 地 拟 合 1:1 线 ， 
锡 林 格 勒 站 点 拟 合 曲线 斜率 达 1.0031 (R°=0.81, P< 
0.05) ,多 伦 站 点 拟 合 曲线 斜率 达 1.0002(CR=0.90,P< 
0.05) ,两 个 站 点 的 BEPS 模型 模拟 结果 与 观测 数据 


(8) 


(b) 多 伦 


y-1.0002x-0.16 
R^-0.90, P«0.05 


20 40 60 80 100 
观测 ET/mm 


图 2 站 点 观测 数据 与 BEPS 模拟 蒸 散 (CE7) 比 较 
Fig. 2 Comparison between the observation data and evapotranspiration (ET) modeled by BEPS 
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基本 吻合 ,表明 BEPS 模 型 能 够 较 好 地 模拟 内 蒙古 
半 干 旱 区 的 ET。 

2.2 蒸 散 (ET) 的 时 空 变 化 格局 

2.2.1 KEET) ARRA ”根据 STARS 突变 检 
验方 法 ,内 蒙古 半 干 旱地 区 ET Te 1981—2018 4E f£ 
在 2 次 突变 ,分 别 发 生 在 1997 年 和 2013 年 ,对 应 RSI 
值 为 1.04.0.49 ,分 别 达 到 了 0.01,0.05 显著 性 水 平 
(图 3a)。CA 曲线 判断 突变 年 份 为 1997 年 ,内 蒙古 
半 干 旱 区 57 总 体 呈 现 * 降 低 - 增 高 "趋势 ,1997 年 前 
距 平 以 负 值 为 主 (-6.29~-37.56 mm-a) ,1997 年 后 
ET 距 平 以 正 值 为 主 (-10.58~52.95 mm - a`') (图 
3b)。21 世纪 以 后 累积 距 平 曲线 的 拐点 并 不 明显 ， 
且 由 于 2013 一 2018 年 跨度 短 , 综 合 两 种 突变 检验 方 
i nf ELA y , 1981—2018 Æ ET T£ 1997 4E Az ^E HH T 
突变 ,因此 ,将 研究 时 期 分 为 1981 一 1996 年 和 
1998 一 2018 年 两 个 时 期 。 

2.2.2 蒸 散 (PET) 的 时 空 分 布 特征 1981—2018 年 
内 蒙古 半 干 旱地 区 ET 以 1.75 mm- a^ I i EE S8 RUE 
动 上 升 趋 热 (R=0.73, P«0.05) ; ET 年 总 量 最 低 值 
(286.14 mm) 出 现在 1992 年 ,2017 年 最 高 达 376.65 
mm( 图 3c)。1981 一 1996 年 有 7 无 明显 变化 趋势 ( 增 


400 
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—e— ET 
350 


ET/(mm:a) 


y71.75x-3184 
R?—0.73, P«0.05 
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速 为 0.41 mm .a', P>0.05) , 1998—2018 4E ET 以 
1.70 mm- a” 的 速度 增加 (R=0.51,P<0.05)。 

由 内 蒙古 半 干 旱地 区 ET 的 空间 分 布 及 变化 可 
知 (图 4a),1981 一 2018 年 内 蒙古 半 干 旱地 区 的 多 年 
平均 E87 呈 中 间 低 两 边 高 的 分 布 特征 。 这 与 研究 区 
的 土地 利用 类 型 密切 相关 ,多 年 平均 ET 总 量 低 于 
300 mm 和 300~350 mm 的 区 域 占 人 研究 区 总 面积 的 
28.17% 和 50.44% , 主要 分 布 在 内 蒙古 高 平原 南部 草 
原 和 未 利用 地 。 多 年 平均 ET 总 量 在 350~400 mm 的 
区 域 主 要 分 布 在 大 兴安 岭南 部 和 松 辽 平原 ,以 林地 
和 耕地 覆盖 为 主 , 占 人 研究 区 面积 的 15.98% ;多 年 平 
均 T 总 量 高 于 400 mm 的 区 域 主要 分 布 在 大 兴安 岭 
南部 和 内 蒙古 高 平原 南部 的 林地 , 占 研究 区 总 面积 
的 5.41%。1981 一 2018 年 内 蒙古 半 干 旱 区 ETME 
趋势 呈 增加 趋势 ,其 中 99.70% 的 区 域 ET 显 著 增加 
(图 4b)。ET 以 0~1 mm-a” 的 速度 递增 的 区 域 零 星 
分 布 于 研究 区 内 , 仅 占 研究 区 面积 的 $.13% ; ET VÀ 


1-2 mm: a 的 速度 递增 的 区 域 占 研究 区 面积 
69.8196; M ET Ai T 2 mm*a 的 速率 递增 的 区 域 占 
人 研究 区 面积 的 24.75%, 主 要 分 布 于 研究 区 东部 和 内 
蒙古 高 平原 南部 的 林地 和 耕 地 。 


累计 距 平 /(mm-a 1!) 
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Fig.3 Mutation and temporal change of evapotranspiration (ET) 
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110?E 120°E 110°E 
T 
(a) ET/(mm-a!) (b) ET 变化 趋势 
ES <300 E <0 
z| EH 300-350 ZL m o- 
$| 国 350~400 +| 8 12 
E >400 m >2 
号 P»0.05 
80% 80% 
z 60% Z | 60% 
sa 40% * 40% 
20% ] 20% | 
300 km oy | 。 30km (vl -可 
图 4 内 蒙古 半 和 干旱 区 莹 散 (E7) 空 间 分 布 特征 及 趋势 变化 
Fig.4 Spatial change and trend pattern of evapotranspiration (ET) 
表 1 蒸 散 CE7) 及 其 影响 因子 在 不 同时 期 的 变化 特征 
Tab. 1 Variation characteristics of evapotranspiration (ET) and its influencing factors in different periods 
1981—2018 年 1981—1996 4: 1998—2018 年 
Bf " 
B R B R B R 
ETl(mm*a') 1.75 0.73 -0.41 0.05 1.70 0.51” 
VPD/hPa 0.00 0.50 0.00 0.03 0.00 0.18 
LAI 0.01 0.15 0.01 0.23 0.03 0.28 
TA/C 0.05 0.45" 0.06 0.22" 0.01 0.02 
WD/(m*s) -0.01 0.64” -0.03 0.74" —0.00 0.02 
SMV/(m?*m?) —0.00 0.03 —0.00 0.01 —0.00 0.002 
PRE/(mm-:a ) —0.82 0.03 2.65 0.12 1.95 0.04 
SD/h -1.24 0.05 -2.01 0.04 -4.91 0.18 


亏 缺 .风速 日 照 时 数 . 叶 面积 指数 和 体积 土壤 水 分 。 下 同 。 


2.3 突变 前 后 蒸 散 (7T) 影 响 因 子 变化 特征 分 析 
从 内 蒙古 半 干 旱 区 kT 及 其 影响 因子 在 不 同时 
期 的 变化 特征 可 知 ( 表 1), 在 1981 一 2018 年 间 ,VPD 
和 TA 分 别 以 0.002 hPa- a^ 10.05 Ca 的 速度 增加 
(P«0.01) , LAI 每 年 以 0.01 的 速度 递增 (P<0.05); 
WD 以 0.01 ms!'…a! 的 速度 递减 (P<0.05)。 在 
1981 一 1996 年 期 间 ,LAIT 和 TA 每 年 以 0.01 和 0.06 *C 
的 速度 递增 (P<0.05);WD 以 0.03 ms1.a! 的 速度 递 
减 (P<0.01)。 在 1998 一 2018 年 期 间 ,VPD 和 LAI 增 
加 速率 分 别 为 0.002 hPa- a^ 10.01 a” (P«0.05) ; SD 
以 4.91 h-a” 的 速度 递减 (P<0.05)。 突 变 前 后 影响 
因子 变化 特征 对 比分 析 , 在 1997 年 (突变 年 份 ) 后 
VPD 表现 出 显著 增加 趋势 ,LAI 增 速 为 1997 年 前 的 
1.90 倍 ,无 论 是 植被 绿化 趋势 还 是 大 气 水 分 亏 缺 的 加 
剧 , 均 有 利于 Z7 增 加 ;但 在 1997 年 后 TA 增加 放 组 ， 


注 :*# ,# 表 示 通 过 0.01.0.05 水 平 显 著 性 检验 ;B 为 变化 趋势 , 尼 为 决定 系数 。TA .PRE .VPD .WD .SD LAIFI SMV 分 别 代 表 温 度 .降水 .大气 水 分 


SD 星 显著 减少 趋势 ,TA 和 SD 表征 了 ET 能 量 供 应 的 
情况 ,其 变化 不 利于 Ek7 增 加 。 综 合 来 看 ,VPD 和 LAI 
导致 的 正 向 影响 主导 了 1997 年 前 后 E87 的 变化 。 
2.4 ET 变化 影响 因子 分 析 

ET 受 多 种 气候 因子 和 植被 共同 影响 。 区 域 整 
体 上 1981—2018 年 ET 受 VPD、TA、WD、SMV 和 
PRE 气候 因素 及 以 LAI 为 代表 的 植被 动态 变化 共同 
影响 ( 表 2)。 其 中 VPD .LALTA 与 £7 间 存 在 显著 
正 相 关 关 系 (P<0.05) ,WD 、SMV 和 PRE 5j ET [H] f£ 
在 显著 负 相 关 关 系 (P<0.05)。1981 一 1996 年 A7 变 
化 不 显著 ,其 变化 受 VPD 和 PRE 共同 影响 (P< 
0.05) ,其 中 TA 和 SMV 的 影响 虽然 未 达到 显著 性 水 
平 ,但 相关 系数 达到 了 0.41 和 -0.49。1998 一 2018 年 
VPD 增 加 导致 了 水 汽 梯度 增加 ,有 利于 ET 增加 (P< 
0.01) ,LAIT 和 TA 与 £7 间 存 在 显著 正 相 关 关 系 (P< 
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表 2 蒸 散 (E7) 与 影响 因子 的 相关 性 
Tab.2 Correlation between evapotranspiration (ET) and 


climatic factors 


PDC BBURHAE EH DLP 1619 


因子 1981 一 2018 年 1981 一 1996 年 1998 一 2018 年 
VPD 0.85” 0.64” 0.70” 
LAI 0.81” 0.11 0.54 

TA 0.70" 0.41 0.45 
WD -0.57" 0.21 0.21 
SMV -0.44^ —0.49 -0.60" 
PRE -0.41^ -0.51* 0.02 

SD 0.08 0.43 0.02 


注 :*** 表 示 通 过 0.01 .0.05 水 平 显著 性 检验 。 


图 5 1981—20184 


FRET) 2E UO TS PT 


Fig.5 Path analysis of evapotranspiration 
(ET) during 1981-2018 


0.01) ;虽然 土壤 水 分 亏 缺 的 加 剧 则 不 利于 ZT 的 增 
综合 影响 导致 了 7 急剧 增加 。 
由 于 相关 分 析 无 法 分 离 自 变量 间 的 相关 性 ,以 
研究 区 ET 数 据 为 因 变量 ,同步 观测 数据 TA PRE, 
WD .SD,VPD,LAI,SMV Jy B 25d 


加 (P<0.01) ,多 因子 


7 


Eri fe 4 pr CER 


5, 表 3)。 考 虑 各 因子 直接 和 间接 作用 的 共同 贡献 ， 
对 ET 变化 的 影响 依次 为 VPD>TA>WD>PRE>LAI> 
SMV>SD ,其 中 直接 作用 依次 为 VPD>TA>WD>LAI> 
PRE>SD>SMV ,总 间接 作用 依次 为 TA>SMV>PRE> 
VPD>WD>SD>LAI。VPD 的 决定 系数 最 大 为 0.95， 
这 与 相关 分 析 的 结果 相同 ,其 中 直接 作用 系数 为 
1.12 ,间接 作用 系数 为 -0.27, 表 明 VPD 对 ET 的 影响 
以 直接 影响 为 主 ;TA 决定 系数 仅 次 于 VPD(-0.42)， 
其 中 间接 作用 远大 于 直接 作用 (1.30>|-0.60| ) ,由 
于 TA 是 导致 大 气 干旱 的 重要 原因 之 一 ,VPD 和 TA 
为 彼此 贡献 了 主要 的 间接 影响 。 而 SMYV 的 间接 影 
响 大 于 直接 影响 (|-0.39|>|-0.05| ) ,主要 来 自 于 TA 
和 VPD。 因 此 ,结合 因子 表征 的 物理 意义 、 相关 和 
通 径 分 析 结 果 ,选择 VPD LAIL TA, WD FI SMV 作为 
817 的 主要 影响 因子 进行 贡献 率 研 究 。 

基于 逐步 多 元 回归 分 析 量 化 了 1981 一 2018 年 、 
1981 一 1996 年 和 1998 一 2018 年 3 个 时 期 研究 区 ET 
变化 的 驱动 因素 及 其 相对 贡献 (图 6)。 在 1981 一 
2018 年 整体 来 看 ,VPD 4i 93.5696 EX EX, ET ZHE 
导 因 子 ,平均 能 够 解释 ET 变化 的 47.23%(24.83%~ 
90.4696 ) ; LAI, SMV, TA 和 WD 分 别 在 3.85% 、 
1.5696 .0.89% 和 0.14% 的 区 域 对 有 7 变化 做 出 主要 贡 
献 。 从 两 个 不 同时 期 来 看 ,VPD TA 和 LAI 对 ZE7 变 
化 的 贡献 较为 明显 , 且 VPD 和 LAI 对 ET 变化 的 贡献 
增加 。1981 一 1996 年 VPD、LAI 和 TA TE 65.7996 , 
15.44% 和 12.74% 的 区 域 主导 了 7T 的 变化 ,平均 能 
够 解释 ET 变化 的 43.61% 、37.84% 和 40.79% ; SMV 
和 WD 分 别 在 4.14% 和 1.89% 的 区 域 对 ET 变化 做 出 
主要 贡献 。1998 一 2018 年 显著 增加 的 VPD 和 LAI 
在 研究 区 内 76.46% 和 17.31% 的 区 域 主导 ET 的 变 


RI 蒸 散 CE7) 影 响 因子 的 通 径 分 析 表 


Tab.3 Path analysis statistics of influencing factors of evapotranspiration (ET) 


. 间接 作用 
因子 TREH 总 间接 作用 决定 系数 
PRE WD SD VPD LAI SMV 
TA -0.60 1.30 - 0.03 0.22 0.01 0.99 0.03 0.02 -042 
PRE -0.11 -0.30 0.15 - -0.03 005 . -053 0.09 | -0.02 0.05 
WD -0.34 -0.23 040 -0.01 - -0.02 | -056 -0.04 0.00 0.19 
SD -0.09 0.17 0.05 007 -0.07 - 0.15 -0.05 0.02 -0.01 
VPD 1.12 -0.27 -0.54 0.05 017 -0.01 - 0.02 0.03 0.95 
LAI 0.17 0.08 -0.11 -0.06 0.09 0.02 0.15 - -0.01 0.04 
SMV -0.05 -0.39 030 | -005 . -0.01 0.00  -0.69 0.03 - 0.02 
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igi bid ECL a TAL H TA E VPD 化 ,平均 能 够 解释 ET 变化 的 45.26% 和 39.63%;TA、 
AN SMV 和 WD 分 别 在 2.14%、2.99% 和 1.10% 的 区 域 对 
80 ET 变化 做 出 主要 贡献 。 
2.5 不 同 土地 利用 类 型 影响 因子 对 ET 变化 贡献 率 
va 对 比分 析 
S 40 基于 逐步 多 元 回归 分 析 量 化 了 1981—2018 年 
5 不 同 土地 利用 类 型 影响 因子 对 B7 变 化 的 贡献 ( 表 
E 4)。 在 5 种 土地 利用 类 型 上 ,VPD 是 BT 变化 的 主要 
0 影响 因子 ,其 次 为 LAI。VPD 在 裸 地 和 草地 主导 ET 
1981—2018 rw 1998—2018 恋 化 的 区 域 最 大 , 分 别 m 其 面积 的 96.14% 和 


注 :LALVPD TA, WD SMV 分 别 代表 叶 面 积 指数 .大气 水 分 亏 缺 、 
温度 、 风 速 .土壤 湿度 。 

图 6 1981—20184F ,1981—1996 ^F ftl 1998—2018 4E ZE TI 
(BE7) 驱 动因 子 占 比 

Fig. 6 The proportions of evapotranspiration (£T) driving 

factors from 1981-2018, 1981-1996 and 1998-2018 


95.1896 ,耕地 占 比 最 小 (86.10% ) ,在 不 同 土地 利用 类 
型 上 平均 贡献 率 均 在 45% 以 上 。LAI 在 耕地 和 城市 
贡献 率 占 比 相 对 较 大 (9.82% 和 8.26% ) ,在 草地 和 耕 
地 占 比较 低 , 仅 有 2.43% 和 1.83% ,在 不 同 土地 利用 类 
型 上 平均 贡献 率 均 在 37.65%~38.61%。 其 他 因子 在 
ET 变化 中 贡献 起 主导 作用 的 占 比较 小 ( 均 低 于 3%)。 


表 4 不 同 土地 利用 类 型 影响 因子 对 蒸 散 (B7) 变 化 贡献 率 分 析 


Tab.4 Contribution rate analysis of influencing factors on evapotranspiration 


(ET) variation on different underlying surfaces 


HF TA/% LAI/% VPD/% WD/% SMV/% 
类 型 占 比 0 占 比 0 占 比 0 占 比 0 占 比 0 
耕地 1.30 40.47 9.82 38.10 86.10 45.83 0.22 35.58 2.56 34.73 
林地 0.69 50.23 4.39 38.61 93.72 45.00 0.26 33.56 0.94 34.19 
草地 0.76 43.89 2.43 37.65 95.18 47.86 0.10 39.29 1.53 37.24 
城市 0.96 38.12 8.26 38.12 88.89 48.72 0.14 35.09 1.75 36.66 
裸 地 1.13 40.64 1.83 37.80 96.14 46.89 0.12 41.99 0.79 40.92 


注 : 9 代表 对 应 因子 的 相对 贡献 率 。 


3 讨论 


3.1 模型 不 确定 性 分 析 

研究 通过 站 点 通 量 观测 Zi7 数 据 对 BEPS 模 型 模 
拟 ET 结 果 对 比 验证 发 现 ,模拟 结果 基本 符合 1:1 线 ， 
决定 系数 在 0.81~0.90。 现 有 研究 表明 ”下 ,基于 站 
点 和 文献 搜集 的 实测 数据 对 多 种 87 模型 结果 验证 
的 决定 系数 多 介 于 0.5~0.8;BEPS 模型 具有 精确 的 
ET 模拟 精度 。1981 一 2018 年 BEPS 模 拟 的 多 年 平均 
ET 结果 介 于 286.14~376.65 mm( 图 3); 牛 忠 恩 等 
基于 PT-JPL 模 型 模拟 了 2000 一 2015 年 中 国 &7 情 
况 ,其 中 内 蒙 地 区 多 年 平均 BT 介 于 300~400 mm; E 
思 如 等 基于 WaVEM 模 型 模拟 了 1964 一 2012 年 科 
尔 沁 沙 地 ET 变化 (200~450 mm) ,以 上 均 与 本 研究 


模拟 结果 基本 吻合 。 张 巧 风 等 5 基于 MOD16 数 据 
反 演 了 锡 林 格 勒 草原 2000 一 2014 年 多 年 平均 E87 为 
218.08 mm, 该 值 低 于 BEPS 模 型 模拟 的 草地 多 年 平 
F ET(328.02 mm) ; 而 有 观测 指出 内 蒙古 典型 草原 
生长 季 ET 介 于 210.40~378.80 mm ,退化 或 荒漠 草原 
的 生长 季 ET 介 于 150.90~207.50 mm ? ,MOD16 产 
品 的 反 演 存在 偏 小 问题 申 。 因 此 ,从 站 点 、 文 献 和 
卫星 数据 结果 对 比 可 以 判断 ,BEPS 模型 模拟 结果 与 
实际 吻合 ,具有 准确 陆地 生态 系统 ET 模拟 能 
3.2 格局 突变 原因 分 析 

为 了 确定 1981—2018 年 研究 区 ET 变化 特征 ， 
研究 经 多 种 突变 检验 方法 共同 确定 了 ET 在 1997 年 
发 生 突变 ,1997 年 后 BT 增加 显著 。jJung 等 ' 在 全 
球 尺度 研究 了 1982 一 2008 年 ET 变化 情况 并 将 1997 
年 视 为 研究 时 期 中 的 节点 ,这 与 本 人 研究 相符 。Jung 
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等 引 将 ET 变化 原因 归结 为 1998 年 发 生 的 厄尔尼诺 
现象 ,而 在 内 蒙古 地 区 表现 为 1998 年 的 特大 洪涝 灾 
害 汪 ,但 这 种 降水 极端 事件 导致 研究 区 日 照 缺 乏 和 
VPD 增 加 ,进而 致使 1998 年 ET7 存 在 一 定 程度 的 减 
小 , 距 平 为 -9.39 mma'( 图 3b, 图 3c)。 同 时 ,也 有 研 
究 者 将 A7 变 化 归结 为 VPD 的 显著 增加 。Yuan 580" 
提出 全 球 海洋 ZT 表现 出 1998 年 前 后 从 上 升 趋势 到 
下 降 趋势 的 明显 过 渡 特 征 ;饱和 水 汽 压 随 气候 变 暧 
增强 ,同时 海洋 输送 到 大 气 的 水 汽 减 少 导 致 实际 水 
汽 压 降低 ,进而 导致 了 VPD 的 显著 增加 ,这 种 现象 
与 本 研究 中 的 VPD 在 1997 年 存在 突变 并 在 1997 年 
后 迅速 增加 的 变化 情况 相 吻 合 ( 表 2)。 同 时 ,中 国 
自 1998 年 包括 研究 区 在 内 的 区 域 实施 了 三 北 防护 
林 体 系 建 设 工程 .退耕 还 林 还 草 工程 等 多 项 生态 修 
复工 程 ,中 国 干旱 半 干旱 区 植被 履 盖 条 件 得 到 了 明 
显 改善 趾 , 研 究 区 内 1998 一 2018 年 LAI 的 迅速 增加 
也 验证 了 这 一 点 ( 表 1)。 

3.3 ET 影响 因素 分 析 

陆地 生态 系统 B87 年 际 变 化 主要 受 能 量 .动力 、 
水 分 供应 大气 水 分 亏 缺 情况 和 生态 系统 状况 等 因 
素 的 有 影响。1981 一 2018 年 期 间 ,E7 在 1997 年 前 受 
VPD LAITA 和 WD 的 共同 影响 ,由 于 各 因子 间 存 
在 相悖 的 变化 趋势 ,7 在 多 因素 的 相互 影响 下 变化 
特征 不 显著 ;而 在 1997 年 后 VPD 和 LAI 的 显著 增加 
弥补 了 TA 和 SD 等 其 他 因素 对 ET 的 不 利 影响 ,ET 
在 1997 年 后 呈现 显著 增加 趋势 (图 3 )。 

VPD 作为 大 气 水 分 气 缺 情况 的 表征 因子 ,直接 
影响 裸 地 和 植被 表面 的 水 汽 梯 度 ,进而 影响 ET, 其 
影响 与 温度 也 密切 相关 。 在 水 分 限制 区 域 ,往往 大 
气 干旱 和 土壤 干旱 同时 发 生 “ ,而 在 内 蒙古 半 干 旱 
区 域 87 变 化 由 大 气 干旱 主导 (图 6)。 持 续 升 温 加 
剧 了 VPD 对 ET 的 影响 中 ,大 气 干旱 加 剧 导致 水 汽 
梯度 加 剧 ,增加 了 植被 和 裸 土 上 的 水 分 流失 ,导致 
区 域 干旱 加 剧 ; 但 由 于 内 蒙古 土壤 水 含量 呈 降 低 趋 
势 ,水 分 散失 将 受到 土壤 干旱 的 制约 纹 。 但 是 高 
VPD 可 能 会 触发 植被 的 自我 调节 机 制 ,导致 气孔 闭 
合 以 减少 水 分 散失 ,在 一 定 程度 上 缓解 大 气 干旱 导 
致 的 ET 骤 增 中 。 这 点 在 研究 中 也 有 所 体现 ,在 耕地 
和 林地 由 于 其 自身 的 调节 作用 更 强 , 对 VPD 的 响应 
小 于 草地 和 裸 地 ( 表 4)。 

近 几 十 年 来 全 球 植被 呈 绿 化 趋势 ,植被 的 生物 
物理 反馈 将 影响 全 球 通 量 分 配 呈 。 植 被 对 7 的 影 
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响 机 制 十 分 复杂 ,植被 通过 蒸腾 作用 为 自身 水 分 吸 
收 和 运输 提供 动力 ,植被 增加 将 导致 落 腾 作用 的 增 
加 ;但 同时 植被 绿化 造成 地 表 谈 挡 和 地 表 粗 烟 度 增 
加 ,不 利于 土壤 蒸发 学 。Jiang 等 "基于 Shuttleworth- 
Wallace 模 型 的 研究 指出 1981 一 2010 年 黄土 高 原 干 
旱 半 干旱 区 蒸腾 能 够 解释 大 部 分 ET 增加 ,蒸发 无 明 
显 变化 ; 且 以 植被 增加 为 主导 区 域 占 E7 和 蒸腾 增加 
区 域 的 63.82% 和 78.47%。 在 内 蒙古 区 域 与 黄土 高 
原 区 域 类 似 , 经 历 了 大 规模 政策 主导 的 生态 工程 ， 
且 研 究 指出 在 内 蒙古 半 干 旱 区 植被 增加 也 是 导致 
ET 增加 的 主要 原因 之 一 ( 表 2)。 在 水 分 限制 地 区 的 
生态 恢复 过 程 中 ,人 工 林 的 增加 加 速 了 地 表 的 水 分 
散失 ,破坏 了 区 域 水 平衡 中 ;因此 生态 工程 应 在 当 
地 水 资源 容量 范围 内 进行 以 使 生态 恢复 更 具 可 持 
续 性 。 


4 结论 


本 文 基于 BEPS 模型 ,利用 遥感 数据 及 气象 观 
测 数 据 ,模拟 了 1981 一 2018 年 内 蒙古 半 干 旱 区 生态 
系统 ET, 并 进一步 分 析 了 ET 的 时 空 变化 特征 及 其 
影响 因素 。 主 要 结论 如 下 : 

(1) BEPS 模 型 模拟 结果 与 基于 通 量 观 测 的 实 
W ET 数据 能 够 较 好 符合 1:1 线 , 尼 在 0.81~0.90 ,能 
够 较 好 地 模拟 内 蒙古 半 干 旱 区 的 ET 情况 。 

(2) 1981 一 2018 年 内 蒙古 半 干 旱地 区 ET REI 
以 1.75 mm: a 3€ RE EJ] E FHEA (P&0.05) ;其 变 
化 存在 具有 年 代 际 差异 ,在 1997 年 发 后 明显 突变 。 
1981—1996 年 87 无 明显 的 增加 趋势 ,1998 一 2018 
4E. ET VA 1.70 mm-a 的 速度 增加 (P<0.05)。 

(3) VPD 是 1981 一 2018 年 有 7 变化 的 主要 影响 
因素 ,在 93.56% 的 面积 上 主导 了 区 域 £7 的 变化 ,能 
够 解释 ET 变化 的 24.83%~90.46%; 通 径 分 析 的 决定 
系数 达 0.95, 以 直接 影响 为 主 。VPD 和 LAI 是 1997 
年 前 后 87 不 同 变 化 趋势 的 主要 影响 因素 ,其 显著 增 
加 弥补 了 其 他 因子 的 不 利 影响 ,导致 了 ET 的 显著 
增加 。 

(4) 在 耕地 、 林 地 .草地 ,城市 和 裸 地 5 种 土地 
利用 类 型 上 ,VPD 均 为 主要 影响 因素 , 主导 面积 在 
86.10%~96.14% ,平均 贡献 率 均 在 45% 以 上 ;LAI 次 
之 ,主导 面积 在 1.83% ~9.82% ,平均 贡献 率 均 在 
37.65%~38.61%。 同 时 ,耕地 和 林地 由 于 其 自身 的 
调节 作用 受 VPD 的 影响 小 于 草地 和 裸 地 。 
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Evapotranspiration changes and its attribution in semi-arid 
regions of Inner Mongolia 


ZHAO Xiaohan", ZHANG Fangmin", HAN Dianchen", WENG Shengheng" 
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Information Science & Technology, Nanjing 210044, Jiangsu, China; 2. Jiangsu Key Laboratory of Agricultural 
Meteorology, Nanjing University of Information Science & Technology, Nanjing 210044, Jiangsu, China) 


Abstract: Understanding the variation characteristics of evapotranspiration (ET) and its influencing factors are 
essential for regional water resources. Based on the boreal ecosystem productivity simulator model, remote 
sensing data, and meteorological data, the characteristics of ET change in semiarid regions of Inner Mongolia 
from 1981 to 2018 were simulated and its driving factors were quantified. ET of semiarid regions in Inner 
Mongolia showed a fluctuating upward trend with a rate at 1.75 mm-a ' (P«0.05) from 1981 to 2018. Further, 
there were concurrent differences in ET with a significant mutation happened in 1997. ET increased at a rate of 
1.70 mm.a in 1998-2018 (P«0.05). Vapor pressure deficit (VPD) and LAI were the main driving factors of ET 
changes after 1997, which significantly increased at a rate of 0.002 hPa and 0.01 per year (P«0.05). The adverse 
effects of other factors were suppressed in positive influence of VPD and LAI, which led to a significant increase 
in ET. VPD was the main driving factor of ET change, which dominated regional ET change in 93.56% of the 
area and explained 24.83%-90.46% of ET change, where the coefficient of determination for path analysis was 
0.95. VPD was the primary factor driving the five land use types of cultivated land, forestland, grassland, urban 
land, and bare land with average contribution rates of more than 45%. 


Keywords: evapotranspiration; vapor pressure deficit; BEPS model; leaf area index; Inner Mongoia 


